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Single integro-differential wave equation for Levy walk 
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The integro-differential wave equation for the probability density function for a classical one¬ 
dimensional Levy walk with continuous sample paths has been derived. This equation involves a 
classical wave operator together with memory integrals describing the spatio-temporal coupling of 
the Levy walk. It is valid for any running time PDF of a walker and it does not involve any long¬ 
time large-scale approximations. It generalizes the well-known telegraph equation obtained from the 
persistent random walk. Several non-Markovian cases are considered when the particle’s velocity 
alternates at the gamma and power-law distributed random times. 

PACS numbers: 


Levy walk is a fundamental notion in physics and biol¬ 
ogy [l| with numerous applications including the trans¬ 
port of light governed by Levy statistics 0, anomalous 
superdiffusion of cold atoms in optical lattices i, T- 
cell motility in the brain [^, endosomal active transport 
within living cells In the last few years, the interest 
in the Levy walk models increases rapidly in the context 
of Levy flight foraging hypothesis according to which a 
Levy transport is optimal to search for randomly located 
objects Levy movement pattern has been observed in 
microorganisms, insects, molluscs, birds, etc. Q- Recent 
publication Levy walk in Review of Modern Physics pro¬ 
vides a detailed discussion of existing applications and 
the most current status of Levy walk theory Q . 

The derivation of the governing equations for a Levy 
walk from an underlying stochastic movement is a long¬ 
standing problem. Many contributions on this subject 
have been given since the pion eering works three decades 
ago (see, for example, [3, [l3])- The standard model for 
a Levy walk is based on the continues time random walk 
(CTRW) with coupled probability density function for 
the jump length and the waiting time between two suc¬ 
cessive jumps. Two integral equations for the probability 
density function (PDF) p{x,t) for a walker position x at 
time t and arrival rate j {x, t) can be formulated and 
solved by the Fourier-Laplace transform technique di- 
An equivalent single integral equation for p(x, t) has been 
also formulated . To describe a Levy walk spatiotem- 
poral coupling, two dynamical equations involving frac¬ 
tional material derivatives have been suggested in [ij . 
Another approach to describe the superdiffusive behav¬ 
ior is based on the analysis of joint PDF p{x, v, t) of the 
particle’s position x and velocity v. Various fractional 
generalizations of Kramers-Fokker-Planck equation for 
p{x,v,t) have been derived [l3l - [l5j . However, to the 
author’s knowledge, the single integro-differential equa¬ 
tion for the Levy walk PDF p{x,t) is still not available 
in the literature. In this paper we obtain the following 
integro-differential wave equation 


where u is a constant speed of walker, ip(u) is the velocity 
jump density: 

t{u) = ^S{u-v) + ^S{u + v) (2) 

in the velocity space V and K (r) is the standard memory 
kernel from the CTRW theory [I|- R is determined by 
its Laplace transform K{s) = 'ijj{s)/^{s), where 'if{s) and 
111 ( 5 ) are the Laplace transforms of the running time den¬ 
sity s-nd the survival function T'(r). From Eq. ([1]) 
with p{x, 0) = po (x) and pt(x, 0) = 0, one can obtain the 
well-known expression for the Fourier-Laplace transform 
of the PDF p(x, t) 


p{k,s) = 


(s + ikv) +(s — ikv) po{k) 


2 — '!/>(s -|- ikv) — %f{s — ikv) 


(3) 


where po(fc) = fj^po (x) e'-^'^dx 

Derivation. We consider the Levy walk as the ran¬ 
dom particle’s motion with continuos sample paths (no 
jumps) along one-dimensional space. The particle starts 
to move with constant speed v at time t = 0 and af¬ 
ter a random time (running time) it either continues the 
movement in the same direction with curtain probabil¬ 
ity or changes the direction and moves with the same 
constant speed. The random running time is defined 
by the switching rate 7 (t) or the running time PDF 
= 7 (r) exp [— fj 7 (s)(is] . To derive the govern¬ 
ing equation for the Levy walk, we start with Marko¬ 
vian model involving structural densities with the extra 
running time variable r [Hill. We define the struc¬ 
tural PDF’s of walker, n+(a:,t, r), at point x and time 
t that moves in the right direction, (-I-) , with constant 
speed V during time r since the last switching. The prob¬ 
ability density function n-{x,t,T) corresponds to the 
walker that moves in the negative direction, (—). The 
balance equations for both structural PDF’s n+{x,t,T) 
and 7T,_(x, t, t) can be written as 


d^P _ 2 ^ , 
dt^ ^ dx‘^ 


[ [ K{t) (p{u) 
Jo Jv 


f d d\ 


X 


( 1 ) 


dn± 


± V 


dn± 

dx 


dn± 


-7(T)n±, 


(4) 


where the switching rate 7 ( 7 ") depends on the running 
time T. If the walker moves in the positive direction it 


p(x — UT, t — t) dudr = 0, 
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can switch with rate 7 ( 7 ") either to the opposite direc¬ 
tion with the probability a_ or keep the same direction 
with the probability such that a+ -I- a_ = 1. For the 
walker moving in the negative direction corresponding 
characteristics are /3+ and /3_. The well-known velocity 
model and two-state model are just particular cases of 
this general two-state model. For example, the choice 
= /3_ = 0 and a_ = ,5+ = 1 corresponds to the 
two-state model [2l|. The probabilities a± = /3± = 1/2 
correspond to the velocity model [l^. We assume that 
at the initial time t = 0 all walkers start to move with 
zero running time t 

n±(x,0,T) = ^Po(x)S(t). (5) 

Here we consider the symmetrical initial conditions when 
walker start to move on the right with probability 1/2 
and on the left with the same probability. The boundary 
conditions at zero running time can be formulated as 

n±(x,t,0) = a± / 'y{T)n+{x,t,T)dT -\- 

Jo 

P± [ ^{r)n-{x,t,T)dT. ( 6 ) 

Our aim is to derive the master equations for the proba¬ 
bility density functions p+{x,t) and p-(x,t) defined as 

p±{x,t)= / n±{x,t,T)dT. (7) 

By differentiating 0 with respect to time t and using 
the balance equations o we obtain 


dp± 

dt 


, , /■* dn± 

n±{x,t,t)TvJ 

-J - J -f{T)n±{x,t,T)dT. 


This equation can be rewritten as the master equation 


= j±{x,t)-i±{x,t), (8) 

where the rates of switching i±{x,t) are defined as 

i±{x,t)= 'y{T)n±{x,t,T)dT (9) 

Jo 


and the arrival rates j±{x,t) are 


j±{x,t) = n±{x,t,0). (10) 

By using the definitions of switching and arrival rates m 
and cni), Eq. (|6l) can be written in the compact form 


j±{x,t) = a±i+{x,t) -k P±i-{x,t). (11) 


In what follows we consider only the simple case of a sym¬ 
metric Levy walk for which a± = j3± = 1/2. In general, 


the probabilities a± and I3± can be useful to formulate 
the impact of the external force or chemotactic substance. 
Substitution of (HU with a± = P± = 1/2 into the master 
equation ([5]) gives 


9p+ , dp+ 
dt dx 


1 , N 1. / N 


( 12 ) 


where the switching rates i±{x, t) can be found as follows. 
By the method of characteristics, we find from (|3]) 

n±{x,t,T) = n±{x ^ VT,t — T<t, (14) 

where '^(t) is the survival function 

4'(r) = (15) 

Note that at r = t we have a singularity due to the 
initial condition ®. Substitution of da into and 
dZ]) together with the initial condition ([5]) gives 

/■* 1 

i±{x,t) = / j±{xTvT,t-T)ilj{T)dT-po{xTvt)i;{t), 
Jo ^ 


p±{x,t) = J j±{xTvT,t-T)'i> (r) dT + ^po{xTvt)'i> (t). 

Applying the Fourier-Laplace transform to these equa¬ 
tions, we can find expressions ior ipx, t) and i-{x,t) in 
terms of p+{x,t) and p-{x,t) as 

i±{x,t) = / K{t)p±{x ^ VT,t — T)dT. (16) 
Jo 

Note that four integral equations i±(x,t) and p±{x,t) 
above can be written in the standard form of two equa¬ 
tions for the PDF: 


p{x,t)=p+{x,t)+p_{x,t). (17) 


and the arrival rate j = -I- j_ with the jump density 

w{z\t) = ^5{z — vt) -\-^6 {z vt) . For the symmetrical 
case j = -k i- we obtain the well-known equations 


p{x,t) 



j{x — z,t — t)w(z|t) 4' (r) dzdr -k 



z)w{z\t)dz'^/ (t ), 


j{x,t) 



X — Z,t — T)w{z\T)'lp{T)dzdT -k 



z)w{z\t)dzip{t). 


Note that these two equations are the starting point for 
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the various studies of the Levy walk 


Our main purpose is to reduce the system (HD), (nsi) 
with (IT51) to a single governing equation for the PDF 
p{x,t) defined by (IT71) . First we introduce the flux [U 


J{x,t) = vp+{x,t) — vp-{x,t). (18) 


Then by adding (1121) and m we obtain the standard 
conservation equation 


dp 



(19) 


Equation for the flux J can be obtain by multiplication 
of (fT^ and (fT^ by v and subtraction 

f) J C^T) 

— +v^-£ = -v[i+{x,t)-i.{x,t)]. ( 20 ) 

By differentiating m with respect to t and (EUl) with 
respect to x and eliminating d'^J/dtdx, we obtain 


Now we need to express the last term in (I^Tl) in terms of 
p{x, t) alone. From (ITTl) and (fTSl) we find the expressions 
for p+{x, t) and p-{x, t) in terms of the PDF p(x, t) and 
the flux J(x, t) 

= , 22 ) 

Substitution of into (HU) gives 
d 

b+{x,t) -i-{x,t)] = 


lj„ 'Or)|^(x-»r,i-r) + -(x + „r,i 

dp . .dp , ., 

+v—{x — VT. t — T) — v—lx + VT, t — r) ctr. 
dx dx 


t) 

(23) 


By using (l23l) together with (IT^ and dSJ, we obtain a 
single integro-differential equation for the PDF p{x, t) 


d'^p 2 
df^ ^ dx"^ 


1 

2 

1 

2 


f\{T) 

Jo 

f K{r) 
Jo 


f d d\ 
\dt ~'^dx) 
f d d\ 


p(x — VT, t — t) 
p(x + VT, t — t) 


d-T — 


dT. (24) 


This integro-differential wave equation for a Levy walk 
is valid for any running time PDF '4>{t). It general¬ 
izes the well-known telegraph equation obtained from 
the persistent random walk with the constant rate of 
switching 7 . By using the velocity jump density m , we 


rewrite Eq. (j24p in the compact form O- It is instruc¬ 
tive to compare the uncoupled CTRW with Levy walk. 
When jumps and waiting times of the CTRW are inde¬ 
pendent random variables, one can convert the single in¬ 
tegral equation for p{x, t) into a master equation (iL Il6l|. 
For the Levy walk involving spatio-temporal coupling the 
integro-differential wave equation dH plays the same 
role as the master equation for the uncoupled CTRW. 

Let us now consider several examples of the running 
time PDF ip (r). 

Exponential running time density. In the Markovian 
case with the exponential running time PDF 

'<P{t)= ^exp(^-^) (25) 

for which '4>{s) = (I -|- Ts)~^ and K (r) = T~^5 (t), we 
obtain from (1^^ the classical Cattaneo or telegraph equa¬ 
tion [I^ 


d^p 1 dp 2 ^^P 


(26) 


This hyperbolic equation ensures that the density profile 
propagates with finite speed v. 

Gamma PDF g{T,2,\). For the biological applications 
it is important to consider a running time PDF that takes 
the maximum value not at zero time as the exponential 
density m 0- For example, it was found that the run¬ 
ning time density for a single bacterium might deviate 
significantly from exponential approximation (^ . One 
example of such PDF is the gamma density 


i/;(r) =g(r,2,A) = A^Texp(-AT) (27) 


with 'ip{s) = {s + X) ^ and K{s) = A^ (2A-I-s) ^. 
The memory kernel in Eq. (l24l) takes the form 

K (t) = A^ exp (—2Ar). (28) 


The advantage of this exponential memory kernel is that 
one can localize integro-differential equation (l24l) by di¬ 
rect differentiation of (E^ with respect to time twice: 


d^p 

W 


-4A 


d^p 

W 


5A' 


dt^ 




dx^ 


^ d'^p ,.dp „, 2 2 9^P 


(29) 


Note that non-Markovian particle’s movement with ve¬ 
locities alternating at Erlang-distributed and gamma- 
distributed random times was considered in [ 2 ^. 
Next we consider the anomalous case involving walker’s 
velocities alternating at power-law distributed random 
times [T^ [T^. [^. 

Anomalous enhanced transport. We consider two 
anomalous cases: ( 1 ) strong ballistic case for which the 
mean squared displacement: < x^ t^ and ( 2 ) subbal- 
listic superdiffusion with < x^ t^~^, where 1 < /i < 2 
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0. In the ballistic case, we use the survival function By using (I32L we obtain 
'h (r) = E[—{t/tqY] with 0 <u < 1 for which the 

mean running time is divergent In this case (s — ikv)^~^ [(s + ikv) p{k, s) — po{k)] + 


tpis) 


1 

l + irosf’ 


0 < fi <1 


(30) 


(s + ikv)^ ^ [{s — ikv)p{k,s) — po{k)] = 0. (36) 


and the Laplace transform of the memory kernel K (r) is 


k{s) 


s-ip^s) 

1 -^(s) 



In particular, for p = 1/2 and po{k) = I, we factorize 
Eq. (l36l) and find 

(s^ + ^ p(fc, s) — I = 0. 


For this kernel the main equation m can be rewritten 
in the different forms by using material fractional deriva¬ 
tives [13, ■ We write it in the form 

where the operator is defined by its Fourier-Laplace 
transform 

EC {L/c^p} = (s — iku)^~^ [(s -I- iku)p{k, s) — po{k)]. 

,( 32 ) 

In the subballistic superdiffusive case, one can obtain for 
small s 


^{s)^l-Ts + ATs>^, l<p<2 (33) 

for which the first moment T = 'I'(r)dr is hnite and 
the second moment is divergent. Then 

k{s)^^{l + As^-^) 

as s —!> 0. From o we obtain the following equation 

with the diffusion coefficient £> = It is clear from 
m and dSl) that ballistic and subballistic cases are 
fundamentally different. In the strong ballistic case 
(0 < /i < I), the integral term is in the balance with clas¬ 
sical wave equation, while for the subballistic superdiffu- 
sive case (1 < /r < 2), the memory term is in the balance 
with the Cattaneo (telegraph) equation. One can per¬ 
form various asymptotic analysis of (EH) 3(iid. (j34l) , obtain 
the pseudo-differential equations for the walker’s PDF 

a ition and determine the shape of PDF profiles 

As an illustration let us find the long-time asymptotic 
solution to Eq. (I5T|) - In the limit tq 0, the evolution of 
the PDF p (a:, t) is determined by the integral term, while 
the Hrst two wave equation terms in (EH are irrelevant. 
The PDF p {x, t) obeys 


p(u)Ll_ pdu = 0. 


(35) 


By using the inverse Fourier-Laplace transform we obtain 
the well-known self-similar profile [H 

p {x, t) = Tr~^ ^ . 

For the arbitrary p from the interval 0 < p < 1, one can 
find from that the solution to (pp) is the Lamperti 
distribution [8| . It would be interesting to use a new wave 
equation for the non-normalizable density problem for 
superdiffusive anomalous transport [s^. Note that our 
equation can be generalized for the case when the particle 
moves with the random velocity: x±{t) = in + 13 {t), 
where B (t) is the Gaussian white noise. 

In summary, we derived the integro-differential wave 
equation for the probability density function for a po¬ 
sition of a Levy walker with continuous sample paths. 
This equation involves a classical wave operator together 
with memory integrals describing a spatio-temporal cou¬ 
pling of Levy walk. It is valid for any running time PDF 
and it does not involve any long-time large-scale approx¬ 
imations. It generalizes the well-known telegraph equa¬ 
tion obtained from the persistent random walk. For the 
Levy walk the integro-differential wave equation dH plays 
the same role as the master equation for the uncoupled 
CTRW P, [l^ . Our technique may lead to the significant 
advances in the extension of the linear Levy walk models. 
Our approach may be particularly helpful when we deal 
with the Levy particles interactions (36| when it would be 
difficult to take into account nonlinear terms within the 
standard approach involving inte gral equations (see the 
similar problems for subdiffusion [371. l38l|L The descrip¬ 
tion of the Levy walk in terms of the Markovian balance 
equations (El) proves to be very useful for the analysis of 
Levy walk with a random death process [^ . Our the¬ 
ory can be useful to formulate the impact of the external 
force or the chemotactic substance on the random move¬ 
ment of particles with finite velocities. This theory can 
be also useful for the implementation of the non-linear 
reactions and development of the theory of wave propa¬ 
gation in reaction-transport sj^ems involving enhanced 
diffusion and memory effects [ll, . 
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